Testing of spin ordering Hamiltonian with ultracold atoms in optical lattices 
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Laser cooling and trapping are now widely used in atomic physics laboratory. Interestingly, cold 
atoms in optical lattices are now used in advanced research to mimic phenomena in condensed 
matter physics and also as a test laboratory for the models of these phenomena. It follows then 
that it is now possible and necessary to advance the atomic physics laboratory by including the use 
of ultracold atoms in optical lattices for instructional contents of phenomena in condensed matter 
physics. In this paper, we have proposed how to introduce into the atomic physics laboratory the 
study of quantum magnetism with cold atoms in a double well optical lattice. In particular, we 
demonstrates how to compare the theoretical parameters of a spin Hamiltonian model with those 
extracted from spin ordering experiment. 
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I. INTRODUCTION 



Cold atoms in optical lattices is the application of two formerly distinct aspects of physics: quantum gases from 
atomic physics [l[ and laser theory from quantum optics Q. The optical lattices are artificial crystals of light, that 
is, a periodic intensity pattern formed by interference of two or more laser beams. As an insight, a pair of these 
laser beams in opposite directions (that is, two orthogonal standing waves with orthogonal polarization) will give a 
""O ' one-dimensional (ID) lattice, two pairs in two opposite directions can be used to create a 2D lattice and a similar three 
pairs in opposite directions will give a 3D lattice. Atoms can be cooled and trapped in these optical lattices. Thus in 
simple form, an optical lattice looks effectively like an egg carton where the atoms, like eggs, can be be arranged one per 
well to form a crystal of quantum matter Q . Though the cold atoms in optical lattices was initially used to investigate 
quantum behaviour such as Bloch oscillations, Wannier-Stark ladders and tunneling phenomena usually associated 
with crystals in a crystalline solid 0, @ , it is the theoretical proposal 0] and consequent experimental realization Q of 
the superfluid to Mott insulator (SF-MI) transition which is an important phenomenon in condensed matter physics 
that has given rise to the possibility of using it as a test laboratory for phenomena in condensed matter physics. 
The success of the SF-MI transition in turn emanates from the laboratory observation of Bose Einstein condensation 
(BEC). The history of BEC began in 1924 when Satyendra Nath Bose first gave the rules governing the behaviour of 
' photon which is the commonest boson. Excited by this work, Einstein in the same year extended the rules to other 
bosons and thereby gave birth to the Bose- Einstein distribution (BED) While doing this, Einstein found that not 
qq ' only is it possible for two bosons to share the same quantum state at the same time, but that they actually prefer 
doing so. He therefore predicted that when the temperature goes down, almost all the particles in a bosonic system 
would congregate in the ground state even at a finite temperature. It is this physical state that is called Bose-Einstein 
. ! condensation. Thus it has always been considered a consequence of quantum effects from statistical mechanics in 
many textbooks as the phase transition is achieved without interactions [1,0. The Einstein's prediction, however, 
was considered a mathematical artifact for sometime until Fritz London in 1938 while investigating superfluid liquid 
helium, realized that the phase transition could be accounted for in terms of BEC. This analysis, however, suffered a 
major set back because the helium atoms in the liquid interacted quite strongly. This was why scientists had to move 
ahead in search of BEC in less complicated systems that would be close to the free boson gas model. Fortunately, 
the breakthrough came in 1995 when the first BEC was observed in rubidium atoms and this was followed by similar 
observations in some other cold alkali atoms such as those of lithium and sodium (see more details in Ref. ([10] and 
a guideline to the literature of BEC in dilute gases in Ref. [ll[ . 

As stated above, the observation of BEC led to the observation of (SF-MI) transition and thereby open the possibility 
to investigate various phenomena in condensed matter physics by mimicking them with ultracold atoms in optical 
lattices. This possibility has led to a deluge of studies (see [l2| for a recent review) as it brings together atomic 
physicists, quantum opticians and condensed matter physicists. One draw back is that even when there have been 
theoretical papers investigating these phenomena with fermionic cold atoms [lH, cold bosons are used in the actual 
experiments [14j for testing spin ordering. This has been overcome by the recent observation of the MI with fermonic 
atoms [l5] . It follows then that the possibility to use cold atoms in optical lattices as a test laboratory for condensed 
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matter physics is no longer a speculative physics. Rather, it has become an aspect of physics with its own methods and 
approaches. Therefore, it has reach a stage when it should start having some introductory impact on our curriculum, 
possibly as applications of optics, atomic physics and simulation of spin ordering Hamiltonians [l6[ in condensed 
matter physics. The purpose of this current study is to present a pedagogical study of investigating spin ordering in 
an isolated double - well - type potential (simply double well (DW)) which can be adopted for instructional purposes. 
For the DW is one of the simplest experimental set ups of optical lattices to study spin Hamiltonians [l7| . This is 
because the system can be completely controlled and measured in an arbitrary two-spin basis by dynamically changing 
the lattice parameters 13j. On the theoretical side, the DW can be considered as two localized spatial modes separated 
by a barrier and consequently be investigated as a two-mode approximation Q. 



II. THE DOUBLE- WELL SUPERLATTICE 



The DW is a ID optical lattice in which the transverse directions are in strong confinement and thus the motions 
of an atom in these directions are frozen out. To create the DW, we start with a standing wave of period d (long 
lattice) so that the potential seen by the atoms trapped in it is 

V^x) = V 1 cos 2 (nx/d) (1) 

where V% is the lattice depth, which is a key parameter for a special lattice potential. 

Next we superpose a second standing wave with period d/2 and depth V 2 (short lattice) on the first one as in Fig. 
1(a). This will lead to a symmetric double- well superlattice (Fig lc) with a total optical lattice [Til [l8l| 

V(x) = Vi(x) + V 2 cos 2 [2irx/d]. (2) 

The configuration and varying of the parameter space (i.e. various parameters) of a Hamiltonian to be tested in 
this superlattice is achieved by manipulating and controlling the depths of the short and long lattices. For example, 
by increasing the lattice depth of long-lattice V±, we could reach from superfluid to Mott-insulator regime, which 
is convenient for studying the few particles phenomena in a local double-well cell. And the barrier height of the 
double- well is controlled by the lattice depth of short-lattice, V 2 . The effective double- well is reached if V\ > 4V 2 . 
Otherwise the minimal points of the optical lattice are the bottom of the long-lattice. This could be seen clearly from 
Eq. ([2]) after we expand the cos(27ra/d) term, 



V{x) = W 2 
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so that cos 2 z > results in V(x) m i n = V 2 , that is, the minimum is one of the long-lattice with a lift of V 2 . On the other 
hand, if V\ < 4V 2 , an effective double-well is created. The minimum can be found at cos 2 (irx/d) = [1 — V\ / (AV 2 )] / 2 . 
Making a power series expansion around the potential minimum, then a single atom of mass m trapped initially in 
any of the well will freely tunneling back and forth with the oscillation frequency of 



u = 1 (IWj-V?) 

d\j 2mV 2 K ' 

Thus the frequency depends on not only the lattice depths V\ and V 2 , but also on the lattice spacing d. Usually, the 
small lattice spacing d is preferred as it leads to a large frequency though this could also be restricted by changing 
the ratio Vi/4:V 2 . This preference also lead to the use of the recoil energy of the short lattice as the unit of the depths 
of the optical lattice 

where A is the wave length of the short lattice. 

For example, in the experiment [3], the depth of the long lattice is 10E r while the depth of the short lattice is about 
6E r < V 2 < 44E r . This gives the oscillation frequency in a range 



7T 119E r 7T /7719£: r 

For 87 Rb and d = A = 765nm, we ca n get E r = 2.596 x 1CT 30 J, which gives 109kHz < A < 326kHz. Lets define 
the harmonic oscillator length a = ^Jhjmuj, then we can readily get 118. 6nm < a Q < 204. 5nm. Comparing to the 



3 



d ( &I2 




FIG. 1: (a) Two standing waves in opposite directions and with periods d and d/2 resulting in (b) a chain of double wells from 
which we can study (c) an isolated symmetric double well or (d) asymmetric double well 



period of short lattice d/2 = 382. 5nm, it implies the ground state wave function is rather localized, which ensures the 
validity of the two-mode approximation. 

Finally, it is pertinent to describe how to create an asymmetric DW (Fig Id) . The potential bias or the tilt A of the 
double-well is introduced by changing the relative phase of the two potentials (i.e. short and long lattices) and this 
can be realized by applying a magnetic field gradient of B' [14( . Consequently, tuning B' gives the potential difference 
between the two potential minima of the DW. We can realize the adiabatic and diabatic operations on the tilt of the 
DW by controlling the increasing speed of B' [l9[ . 



III. TWO-MODE MODEL WITH ONLY ON-SITE INTERACTION 



Within the above consideration and assuming atoms confined in an isolated DW, we reach the two-mode approxi- 
mation represented as a two-site version of the Hubbard model [201 ] 
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where a^ L r {<1sl,r) is the creation operator (annihilation operator) for an atom with spin ct(ct) =t (|), | (T)> n a.L.R 
is the corresponding number operator, J (both J and t are used in the literature though the cold matter community 
seems to prefer J) describes the tunneling rate between the two wells, A is the potential bias for the double-well and 
U is the two-body interaction when two atoms occupy the same site. 

Basically, the Hubbard model is a single band lattice model supporting a single atomic state which can hold up to two 
particles. If we consider that these particles have two internal spins, cr(a'), then the Hamiltonian will consist of the 
superposition of six Fock basis states [2l], [22j denoted by | f J,, 0), | j, J.), | |, |), |0, t j), | T, T) an d | 1,1), with a basis 
state \l,r) denoting the I = l eft and r = right wells. The use of the superposition principle which is a fundamental 
concept in quantum theory [23| is consistent with both the model Hamiltonian and experiment in which the time 
evolution of the initial states produces coherent superposition of states. Thus the wavefunction of the system will be 
the superposition of all possible states 
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where for convenience, \<j>u)t,<t = l*)/l s )ff,ff = I 11, 0) ± |0, IT)/V2, \<kj 
\^ij)t,<7 = \t)ta = I t>T) ± I l,l)/V2, with i, j denoting the sites while | s 
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\t)/\s)%a - I T,4>±l 4,T>/V2 and 
| t) denote singlet and triplet states 



respectively. From basic physics, the orientation of the two spins in a state can either be singlet | s) if S = or triplet 
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A. eigenenergies and eigenstates at A = 



By tuning the potential bias A = 0, we can obtain all the eigenenergies and corresponding eigenstates analytically. 
This is achieved by directly diagonalizing the Hamiltonian in Eq. ([7]) to obtain eigenenergies and eigenstates as shown 
in Table (p}. In the weak interacting case, U « J, the state | H, 0) and state |0, ||) have lower energy so that the 
doubly occupied singlet state, \4>u)^ a = I s ) w ^ be the ground state of the system. This can be considered as the 
signature for a superfluid state for bosonic atom in a double well. However, the strong interaction regime is more 
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interesting to study for spin ordering. In this regime, U » J, the ground state will be singly occupied as the large 
atomic repulsion energetically suppress the double occupancy. Here it is the that are occupied while the 

\<l>ij)a a are unpopulated OH^]. The populated | s) and | t) of \(j>ij)^ 9 are nearly degenerate because the energy 
difference between them is about AJ 2 /U, which is a small quantity. However, when J«0, the ground state approaches 
| s) while the first excited state is | t). If we prepare the initial state as antiferromagnetic, | f, J.) = (\t) + |s))/v2, the 
dynamical evolution involves two frequencies [141 ] 



V16 J 2 + U*±U ... 
"^1,2 = 2 " ^ ^ 

From the above frequencies, one could get the tunneling rate J = 1 / 2h^/uJiuj2 and the interaction strength U = 
H(u>i ~ UJ2) respectively. These two frequencies can be obtained from exerimental data and then used to test the 
validity of the simple two- mode model. This has been done for bosonic atoms in experiment [l4T |. The extension 
to fermionic atoms may be different but the eigenenergies and corresponding eigenstates are the same, which means 
we can get similar dynamics as long as the interaction between atoms satisfy U > [13]. On the other hand, the 
interaction of fermion could be also attractive generally. It is interesting to identify the ground state in this situation. 
Table Q works here too. When U < 0, the ground state does not change too much at weak tunneling. However, 



TABLE I: Eigenstates and Eigenenergies of Hamiltonian (J7J) 
Eigenenergy 11,1) ITl.O) IT, I) I I, T) |0,tlHTTT 
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the first excited state is not the triplet state \<j>ij)t s — \t) anymore but the state \4>u)t a = \t)- It is interesting 
that the energy difference is the same AE « 4J 2 /\U\ although the interaction is attractive and the first excited 
state is changed. This analysis shows that \4>u)t s — 10 could be involved in the dynamics if we start from the 
antiferromagnetic initial state | |, 1). The two frequencies that can be observed in experiment are 

V16J 2 + U 2 - U 

hw 1 = u.,hhj 2 = (10) 

Here the two-body interaction strength directly relates to to'-y , which can be extracted from the measured experimental 
data. The results show that ultracold atoms trapped in the superlattice not only could be used to simulate the phe- 
nomena in condensed matter physics, but also offer the possibility to compare the results with theoretical calculation 
of model Hamiltonians. 



B. eigenstates and eigenenergies at A 7^ 



One of the advantages of the trapped ultracold atoms is that they can be controlled precisely and easily. By tuning 
the two optical lattices, we can change the bias of the double-well superlattice. The ground state of the two fermions 
are trapped in the same site in the large potential bias. By slowly reducing the bias, the ground state is followed 
adiabatically to the singlet state. In this way, we can prepare the initial states either in the state | T|,0) or |0, fj,). 
So it is interesting to investigate how the bias influences the states. 

When the potential bias A is included, the eigenstates and the eigenenergies have complicated expressions. Three 
of the eigenenergies are always. The others are the roots of the algebra equation 

x 3 - 2Ux 2 + (U 2 - 4J 2 - A 2 )x + \UJ 2 = 0. (11) 

We numerically solve the equation and plot the eigen spectra in Fig. When U ^S> J, the ground state energy is 
negative and modified by the presence of potential bias. When A is not too big, the ground state energy is close to 
the one without potential bias. Also, the energy difference between the ground state and first excited state is small. 
At large potential bias A 3> U, however, the approximate ground state energy reads 



E g = U- ^4 + A 2 . 



(12) 
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FIG. 2: (Color online.) Spectra of the eigenstates at different potential bias A/ J. Left, U/J= 10.0 in the strong interaction 
regime. Middle U/J = 0.1 in the weak interaction regime. Right, attractive interaction U/ J — —10. 



The two atoms are in the right well and the ground state reaches |0, ||). On the contrary, the system is degenerate 
further in the weak interaction regime U <^ J. If the interaction is attractive, the energy spectra is reversed. The 
energy difference between the ground state and the first excited state is bigger except at A = 0, i.e. an anti-crossing 
appears. This is shown obviously in Fig. The observation from this analysis is that the potential bias can be 

used to control the energy difference between the singlet and triplet states [Hj]. This is why an attempt was made in 
[l4| to use it to drive Eq. into superexchange interaction observed experimentally. The outcome, however, is that 
inter- well interactions have to be included to Eq. [7] to get close to the experimental data. Thus in the next section, 
we will consider such an extension. 



IV. LONG RANGE INTERACTIONS AND INTERACTION INDUCED SPIN-FLIP 



It is obvious from the preceding section that a Hamiltonian to study spin ordering in the cold atoms in optical 
lattices needs to contain long range interactions. It is important to point out that the overlapping of different electronic 
orbitals gives rise to the interaction between spins in condensed matter but this overlapping is very small in optical 
lattices [26| . However, the possibility of the atoms to tunnel through the barrier in quantum mechanics enables the 
inter-site interactions [l2j . Two natural candidates are the inter-site Coulombic interaction V and exchange interaction 
J ex . Interestingly, the inclusion of these interactions as means of going beyond the standard Hubbard Hamiltonian 
to account for ferromagnetism in metals have been proposed p?! |28| . Furthermore, we do not need a potential bias 
since the spin ordering is induced by these interactions. Within these considerations, the extended form of Eq. [7] with 
A = is 



H= E [-J{al L a 9R + al R a sL ) + U(n^ L n lL +n 1R n lB ) + V(n^ L n lR + n lL n^ R ) + J ex ^ ^l^sr^R^l 



(13) 



It is then easy to obtain the ground state energy and wavefunction using the highly simplified correlated variational 
approach (HSCVA) in [2!|. The beauty of this pedagogical approach is that the ground state energy clearly depicts 
the physics of the model as one vary the parameter space as in experiments with optical lattices. Interestingly, the 
method allows the decoupling of the kinetic part from the interaction parts so that we can observe the effects of 
including each of them to the kinetic part. Thus the combination of these two factors makes the HSCVA very suitable 
to investigate the spin Hamiltonian to be tested using cold atoms in optical lattices. 
We start with the variational ground state energy 



E 3 = (W) 



where the H is the model Hamiltonian and the ket in the Hilbert space is the trial wave function (cf. Eq. ([8])) defined as 



iv> = + E x Mh+ E Y *\4i&- ( 15 ) 

i—l,r <i,j>—l,r <i,j> — l,r 

The X and Y in Eq. (fT5|) are the variational parameters. It is straightforward to show [2i| that Eq. (fT4|) leads to a 
3x3 blocked matrix of 2 x 2 and lxl resulting in the lowest state energies [27| . E g /t = E s for the singlet states 
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\s) depending on U, 



and Eg/t = E t for the triplet state \(j>ij 




,± - 
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E t =V- J e . 



(16) 



(17) 



The smallest of these two energies will be the ground state energy of the system. The corresponding eigenvectors 
are then substituted as the variational parameters in Eq. (|15l) to give the corresponding ground state wavefunctions. 
Thus when E s < E tl the system will be antiferromagnetic while it will be ferromagnetic otherwise. Taking into 
account this condition and Eqs. (fl6|) and (fl7|) , the critical value of J ex at which there is transition from one state to 
another is 



Jex 

4J 2 




(18) 



Now to test this Hamiltonian in a double well, we need to know how the atomic positions and spin orientations varies 
with the parameter space. For example, as demonstrated in subsection (A), the ground state of the system will be 
a Mott insulator when the U is very strong. This generally accepted property of the half-filled standard Hubbard 
Hamiltonian (i.e. V = J ex = 0) is already achieved with ultracold fermionic atoms [ll|. One of the signatures of the 
MI state is the decrease in doubly occupied states in the ground state as U is increased. This is demonstrated in Fig. 
3 showing the level of occupation of the states denoted by the variational parameters of the ground state wavefunction 
with increase in U. The inclusion of V , however, enhances the double occupancy and is therefore expected to suppress 
the observation of the MI especially for low values of U. It follows then that when we switch on the J ex , the U is 
likely to drive the system into more singly occupied states and thereby enhancing the transition to a ferromagnetic 
state while the V will suppress it. This is demonstrated in Eqs. (fH)|) - (fT5|) and then depicted in Fig. 4 showing the 
variation of the antiferromagnetic- ferromagnetic transition critical point of J ex with U at various values of V . 
The above theoretical V and J ex can also be compared with the ones obtained from extracted data from the exper- 
iments as was done for J and U in the standard Hubbard Hamiltonian. This is by expressing the possible dynamic 
evolution frequencies from for the singlet states and triplet states as 



hwz 4 = 2 J 




J ex 

4J 



, U V J ex 

1 ± h h — 

47 47 47 



(19) 



Hu 6 = V-J ex . (20) 

We see immediately that we can recover Eq. © from Eq. (|19[) when V = J ex = 0. Taking into account Eqs. ([9]), (fl9l 
and (|20[) , we can then estimate V — i[(<^3 —^4) — (u>i — W2) +^5] and Jex = § [(^3 —^4) — (^1 — ^2) — Thus we 
can also obtain the inter-site interaction parameters from the data extracted from the experiments. 



V. CONCLUSION 



The increasing advancement on how to prepare, manipulate and detect phenomena in condensed matter physics 
using cold atoms in optical lattices has reached a stage when it can be used as instructional means. The fact that 
laser cooling and trapping are now widely used in atomic physics laboratory [30j means the realization of the double 
wells experiment can also be achieved. The first investigation is to mimic the Mott insulator state. By extracting u>\ 
and 0J2 from the experiment, the experimental values of J and U can be compared with the ones from their theoretical 
values. The experiment can then be advanced to determine V and J ex and then compare them with their theoretical 
values as depicted in Fig. 4. Interestingly, the model Hamiltonian studied here has been proposed to account for 
spin ordering in transition metals. It is hoped therefore, that the testing of this extended Hubbard model can easily 
be compared to available data for the transition metals [28| after some refining of the approach here. This will also 
include extending the study to dynamic properties of the model. 
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FIG. 3: (Color online) The decrease of the doubly occupied states Xa while there is increase in inter-site states Xij as the 
on-site Coulomb interactions U increases. 
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FIG. 4: (Color online) The variation of the antiferromagnetic-ferromagnetic transition critical point of the exchange interaction 
J ex with the on-site Coulomb interactions U at various values of the inter-site Coulomb interactions V. A similar graph was 
obtained by Ref. ([ItJ) using different analytical method 
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